Minimal Work Principle and its Limits for Classical Systems. 
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The minimal work principle asserts that work done on a thermally isolated equilibrium system, 
is minimal for the slowest (adiabatic) realization of a given process. This principle, one of the 
formulations of the second law, is operationally well-defined for any finite (few particle) Hamiltonian 
system. Within classical Hamiltonian mechanics, we show that the principle is valid for a system of 
which the observable of work is an ergodic function. For non-ergodic systems the principle may or 
may not hold, depending on additional conditions. Examples displaying the limits of the principle 
are presented and their direct experimental realizations are discussed. 
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Thermodynamics originated in the nineteenth century 
as a science of macroscopic machines, constructed for 
transferring, applying or transmitting energy P]. The 
main results of this science are summarized in several for- 
mulations of the second law . Last fifty years witnessed 
progressive miniaturization of the components employed 
in the construction of devices and machines 0. This 
will open the way to new technologies in the fields of 
medicine, computation and renewable energy sources. 

For microscopic machines and devices, we need to un- 
derstand how the second law applies to small systems. 
There are two aspects of this program: i) Emergence of 
the second law, where one studies fluctuations of work 
or entropy, knowing that on the average they satisfy the 
second law. Important contributions to this topic were 
made by Smoluchowski and others, nearly hundred years 
ago [2] • Recently, the activity in this field revived in the 
context of fluctuation theorems ii) Limits of the sec- 
ond law, where the very formulation is studied from first 
principles. Here we continue on that and study (limits 
of) the minimal work principle based on classical me- 
chanics. This formulation of the second law is opera- 
tionally well-defined for finite systems, and it relates the 
energy cost of an operation to its speed. The principle 
was deduced from experience and postulated in thermo- 
dynamics, where it is equivalent to other formulations of 
the second law. Its derivation in statistical physics was 
formulated several times, on various levels of generality 
0, 0, 0| ■ Almost all these studies concentrate on macro- 
scopic systems and confirm the validity of the principle. 
In Ref. the principle was derived for a finite quan- 
tum systems, and limits, related to energy level crossing, 
were indicated. However, the reduction level of quantum 
mechanics is not always needed, e.g., certain aspects of 
nanoscience are adequately understood already with clas- 
sical ideas In addition, it is not easy to design ex- 
periments in the quantum domain that would check the 
validity of formulations of the second law. Thus it is nec- 
essary to understand the principle in classical mechanics 
and this is our present purpose. Our results will apply 



to small systems, having one degree of freedom. 

Consider a classical system with N degrees of free- 
dom and with a Hamiltonian H(q,p, Rt), where q = 
(gi, qisi) and p = (pi, ...,pn) are, respectively, canoni- 
cal coordinates and momenta. The interaction with ex- 
ternal sources of work is realized via a time-dependent 
parameter Rt = rit/Tn) with time-scale r/j. The mo- 
tion starts at the initial time i; — tiTji, and we follow it 
till the final time tf = Utr- We shall denote z = {q,p), 
dz = dqdp and i?k = H{z,R\^) with k = i, f, for the 
initial and final values, respectively. 

The system starts its evolution from an equilibrium 
Gibbs state at temperature T = 1//3 > 0: Vi{z) ^ 
g-/3ff.(z)/2:i, where = j dze"^^''^). The work done 
on the system is the average energy difference Q 

j dz[Vi{z)Hdz) -V,{z)mz)] (1) 



dsRsV{z,s)dRH{z,Rs), (2) 



where Vi^z) is the final distribution, and where the equiv- 
alence between and is established with help of the 
Liouville equation for V and the standard boundary con- 
dition 'P{z) = Q for z ±oo 0. Eq. 10) justifies to call 
w{z,R) = diiH{z,R) "the observable of work". Con- 
sider processes with different tr, but the same r{t) and 
Rk — f(tk) (k = i, f)- The minimal work principle claims 



AW = W -W >0, 



(3) 



where W is the work for the adiabatically slow realiza- 
tion, where tr is much larger than the characteristic time 
Ts of the system. Eq. (PJ is an optimality statement: 
the smallest amount of work to be put into the system 
{W > 0) is the adiabatic one, and the largest amount of 
work to be extracted from the system {W < 0) is again 
the adiabatic one. The condition of being slow is purely 
operational. If there exist several widely different char- 
acteristic times, there are several senses of being slow. 

Our derivation of the principle consists of three steps 
and assumes that w{z,R) is an ergodic observable of the 
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dynamics with R =const. Note that the known inequal- 
ity W > F[- Fi, where Ff,i = -Tin/ dze'^f '^'^)/'^ are 
free energies referring to the initial temperature T, does 
not in general provide any information about the prin- 
ciple, since for finite systems the adiabatic work is not 
equal to F[ — Fi; see 0111 for examples, and also below. 

1. The initial distribution Vi{z) can be generated by 
sampling microcanonical distribution Ai with initial en- 
ergy probability Pi{E), 



M{z,E,Ri) 



S[E-Hiz,Ri)], (4) 



Vi{z) = J dEP,{E)M{z,E,Ri). (5) 

Hamilton's equations of motion imply ^H{zt,Rt) = 
RtdBH{zt,Rt)- On times ^ r ^ r^; we have for 
the energy change 

ArE=-[H{zt+r,Rt+r)'Hizt,Rt)]= (6) 

T 

(z„ i?,) = — / ds —{z„Rt) + o( — ). 



T ds 



dR 



The last integral refers to the frozen-parameter dynamics 
with Rt — R. Now recall the Liouville theorem dz = dzt 
and energy conservation H{zt+r) — H(zt) — Et, so that 



1 



t + T 



dzw{z)Miz,Et) ^ - J ds J dz w{z)Miz, Et){7) 
dztM{zt,Et)- dsw{z[s;zt]),{8) 



where z[s;zt] is the trajectory at time s, that starts at 
zt at time t. If w{z, R) = djiH{z, R) is an ergodic ob- 
servable of the R =const dynamics, then for r 3> T5 
the time-average in ((S)) does not depend on the initial 
condition zt [IJIj so the integration over zt in © is triv- 
ial, and we get from Q that the time-average in © is 
equal to the microcanonical average at the energy Ft: 
ArE = RtJ dzdRH{z,Rt)M{z,Et,Rt). This implies 
conservation of the phase-space volume f2: 

ArVl{E,R)=Q, n{E,R)= J dz6{E ~ H{z,R)). (9) 

Given the initial energy Ei and the initial and final 
values of Rt, Ef is found from ri(i?i,i?i) = Q{E{,Rf). 
Since dE^{E,R) = uj{E,R) > 0, we can define for 
fixed Rf and i?i two monotonous functions: E{ = (f>{{Ei) 
and its inverse Ei — (j)i{Ei). Note that (pi{E) = 
uji{E)/LUi{(j)i{E)) > 0. The conservation of the phase- 
space volume for ergodic systems is known for almost a 
century . It led to the development of microcanonical 
thermodynamics with entropy defined as In fl |lo| . The 
precision of the above conservation was studied in \vi\ . 



2. Let P{E\E') be the conditional probability of hav- 
ing energies E and E' at t = t{ and t = ti, respectively: 

P{E\E')Pi{E') - (10) 

J dzdz'S[E - Hi{z)]S[E' - Hi{z')]V{z\z')Vi{z'), 

where V{z\z') is the phase-space conditional probability 

V{z\z')=U5[p-p']S[q^q'], U = ^^'''^'^'^-'\ (11) 

and where C{t) = dqH{t) dp — dpH{t) dq is the Liouville 
operator, with e being the chronological exponent. We 
see that / dz''P{z\z') = 1, since U is a. sum of 1 and 
differential operators. This implies together with (0 110() 

P{E\E') = Jdzdz'S[E-H,{z)]^-l^^-^^nz\z'), 

j dE'uji{E')P{E\E')^ujf{E), (12) 

in addition to the normahzation J dEP{E\E') = 1. For 
the adiabatic situation we get from conservation of fi: 



P{E\E') = 5[E - (j)i{E')]. 
3. Now recall (QJ and write AW as 



(13) 



AW^ j dzHi{z)[Vi{z) - Vi{z)] = (14) 
j EdE j dE'Pi{E')[P{E\E') ^ P{E\E')] = J dEgs- 

Integrating by parts, denoting ce{E') = duP{u\E'), 
using (1121), and employmg 0[E ~ (l)f{E')] = e[(l)i{E) - E'], 

gE = j du j dE' Pi{E')[P{u\E')- P{u\E')] 



J " dE'PiiE') - J dE'Pi{E')cE{E') = 



'^''^''\E'PiiE')[l-CE{E')]~ f dE'Pi{E')cE{E'] 

J<jji(E) 

We now employ (|3J) and then (|12|l to obtain 

r<t>dE) 

Z-.gE > e-^^-^^^ J dE'LJi{E')[l - ce{E')] 
_e-/30.(£) f dE'uji{E')cE{E') 



r4>i{E) rE 

/ dE'ojiiE')- / dE'uj{{E') 



■ (15) 



Recahing that uj{E) = fl'{E) and that n{E^in) = for 
the lowest energy i^min, we get (fT^ = 0, i.e., gs > 0. This 
means, from (I14|l . that the principle is proven. Note that 
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FIG. 1: Work W versus the time-scale tr for an ensemble of 
harmonic oscillators with Rt — —b + t^/r^ at initial temper- 
ature T = 2. Full line: b — 0, ti = —tf = tr; the adiabatic 
work is W — 4, as given by Eq. 12311 . Dashed line: b = 0.5, 
ti = —Tr and t[ — tr^/OA. Bold line: 6 — 0.5, = —tr 
and = 0. For the last two cases the adiabatic work (the 
limit Tr oo) is, respectively, plus and minus infinity. The 
minimal work principle, W(tr) > W{oo) for all tr, is seen 
hold in the third case, but not in the first and second case. 



the same proof applies for Pi{E)/u!i{E) being a mono- 
tonically decaying function of E. Neither the proof, nor 
the principle itself, applies to the initial microcanonic 
situation, where Pi{E) = S[E - E-,] flst - 

For non-ergodic systems, where under driving the sys- 
tem can move from one ergodic component (of the 
R =const dynamics) to another, the above proof of the 
principle is endangered, since in general f2 in @ is not 
conserved. Indeed, the argument expressed by QIHl) may 
not apply, since now the time-average in ((SJ depends on 
the ergodic component to which the initial condition zt 
belongs, and in general it cannot be substituted by the 
microcanonical average over the full phase-space. How- 
ever, dnHiz, R) may be ergodic, even though the sys- 
tem is not |l2l| . For an example consider the symmetric 
double-well: H = ip^ — Rtx'^ + gx'^, with g > 0. For 
E <Q, there are two ergodic components related by the 
inversion x —x, but dRH{z, R) = —x"^ is degenerate 
with respect to them. Though the Rt =const motion on 
the separatrix E = Q has an infinite period (due to un- 
stable fixed point q = Q), when the initial distribution is 
microcanonical, the fraction of particles trapped by the 
separatrix is neglegible [l^l, so that for the ensemble ts 
remains finite. Thus ^{E, R) (with the integration in 
over the whole phase-space) is conserved 13]. Thus for 
this case the above proof of the principle applies. 

Limits of the principle. The above derivation of the 
principle assumes that the frozen-parameter dynamics 
supports the microcanonical distribution. This need not 
be always the case. Consider the basic model of the para- 
metric oscillator: H = ^p^ + ^RtQ^. If Rt is always posi- 
tive, the phase-space volume is conserved, and the above 
construction applies. But what happens if Rt touches 
zero at one instant? This is another non-ergodic exam- 
ple, since for Rt = the frozen-parameter phase-space 
consists of two ergodic components with, respectively. 



p > and p < 0. Though diiH{z, R) = is degener- 
ate with respect to them, the microcanonical distribution 
does not exist for Rt — 0. 

To separate the effect of initial conditions, we write the 
solution of the equations of motion q + Rtq — as 

p{t) = OppPi + Opq qi, q{t) = 9qq qi + Oqp Pi, (16) 

where 6ki — dki{t). With Gibbsian initial distribution 
Vi{p, q) cx exp[— ^ — "^f^]' the work reads from |^: 

Next we consider an exactly solvable situation Rt = 
/t]^ (see below for generalizations). The equation of 
motion is solved by substitution q(t) — ^/ix{t): 

q{t) = ci/^ J_i/4(t^) + Ji/a{^), (18) 

2Tfl ^\t\ Ztr 

where J±i/4 are the Bcsscl functions, ci and C2 are to be 
found from initial conditions, and where q{t) is written 
in a way that applies to i < 0: noting Jp{x 0+) = 
2i^r{'p.+i) 01 1 see that q{t) ^ Ci + C2t near t — 0. Let 
us define ti = —TR^/Ri and t{ — tr^/Ri. In the slow 
limit Tij > 1 we need J^{x) = y^cos (a; - ^ - f ) + 
0(x-3/2) Q. Eqs. IHIlHl) then imply 

Oqq = rI^'^ R;^^^{-smu + V2cosv), (19) 

Oqp = Rr^^"^ R^^^*{cosu + V2smv), (20) 

Opq = rI^'^ R^^icosu- V2smv), (21) 

Opp = R^^^'^ Ry^{smu + V2cosv), (22) 

where u — \tr{-/Ri — \fRi), v — ^tr^s/R^ + ^/R{). In- 
serting 1^11211) into lfT7|) we get for the adiabatic work 

W ^T{-l + 3y/Rf/Ri). (23) 

Note that if Rt is always positive, the conservation of 

— 1/2 — 

the phase-space volume oc E(t)Rt gives W — 
{El - Ei) = {^Ih/Ri~ l){Ei) = T[^/lh/Ri- 1). Thus 
Q, is not conserved, if i?t = at one instant. Note that 
the adiabatic^imit still exists, since for a large tr, W 
converges to W; see Fig.^ Eq. shows that the min- 
imal work principle does not hold. Indeed, for a sudden 
variation the Hamiltonian changes, while the state of the 
ensemble does not. This brings 

W.^J dzVi{z)[Hdz) Hi{z)] = T - . (24) 

Ws is sometimes smaller than (|23J, e.g., take Ri = Ri; 
see also Fig. ^ A qualitative picture behind this is that 
for Rt = the particle runs to infinity, and to confine it 
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back (for i?f > 0), the work to be spent is larger for the 
slow case, since for a quick process the particle does not 
have time to move very far; see H24II . 

An extended setup Rt = —b + t'^/r^, b > 0, can 
be worked out with help of hypergeometric functions. 
Fig. n illustrates that the principle is satisfied if Rt de- 
cays monotonically: ti < tf < 0. It is violated if the 
change of Rt is non-monotonic: ti < < t{. The work 
does not saturate for th ^ oo if Rt becomes strictly 
negative. We checked numerically that i) \i Rt is posi- 
tive, but becomes very small, R — 5, the system has two 
widely different characteristic times ~ 1 and Tg that 
diverges for (5^0. Thus there are two senses of being 
slow: the principle is limited for <C t^j <C Tg, while 
it is satisfied if Tg <C t^, since then the phase-space vol- 
ume is conserved, as discussed above, ii) The limits of 
the principle exist for any potential U{q, Rt) that is glob- 
ally confining [U (q) — > oo for q — )■ ±oo) but looses this 
feature for one value of Rt. These limits obviously exist 
for uncoupled particles. We expect that they extend to 
coupled particles put in a (de)confining potential. 

Note that whenever the principle gets limited via the 
above scenario, the slowest process is irreversible. Recall 
that a process is reversible, if after supplementing it with 
its mirror reflection (the same process moved backwards 
with the same speed), the work done for the total cyclic 
process is zero Jj. As seen from (|23|) . the work (equal to 
2T) does not vanish for the cyclic adiabatic process with 
Rt touching zero. Thus the process is irreversible. This 
fact contrasts the quantum limits of the minimal work 
principle found in 13- Those limits are related to energy 
level-crossing, where the adiabatic work is reversible [3- 

Here is an experimentally realizable example that can 
demonstrate the above limits. The simplest LC circuit 
consists of capacitance C and inductance L (the resis- 
tance is either small or compensated) The Hamil- 

tonian is iJ = fj; + where Q (coordinate) is the 
charge, and where $ = (momentum) is the mag- 

netic flux. The parametric oscillator with i?^ — > corre- 
sponds to a time-dependent C (or L) becoming very large 
at some time. Rt becoming negative at some time, can 
be achieved via a negative capacitance Ca < given by a 
special active circuit |l7l | . If such a capacitance is sequen- 
tially added to a positive capacitance C„, then the result- 
ing inverse capacitance = + can be made 
zero and then negative by tuning Cp. The same effect is 
obtained via a negative inductance [ig added in parallel 
to a normal one. As the negative inductance/capacitance 
emulators are widely appplied in compensation of para- 
sitic processes and for improving the radiation pattern in 
antennas, 0, 0, 0| they can serve to test our predic- 
tions. 

In conclusion, we studied the second law in its mini- 
mal work formulation for classical Hamiltonian systems. 
It was shown to hold under the assumption that the ob- 



servable of work (i.e., the derivative of the Hamiltonian 
w.r.t. the driven parameter) is an ergodic function. The 
result applies to small systems. There are, however, nu- 
merous examples of non-ergodicity both for finite and 
macroscopic systems. For such systems we explored sev- 
eral possiblities met in the one degree of freedom situa- 
tion. The minimal work principle applies if the observ- 
able of work is degenerate over ergodic components and 
if the microcanonic equilibrium exists for all values of the 
driven parameters. If the latter condition is not met, the 
principle gets limits. The simplest example of the latter 
is provided by a parametrically driven harmonic oscil- 
lator whose frequency passes through zero. As we saw, 
this situation can be realized experimentally in LC elec- 
trical circuits. Multi-dimensional systems provide more 
complex examples of non-ergodicity. The understanding 
of the second law for such systems still deserves to be 
deepened, in view of the importance of non-ergodicity in 
processes of measurement and information storage [l2j . 
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